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Making use of the Rugh's micro-canonical approach to temperature we study X F-spin systems 
by the strictly energy conserving over-relaxation algorithm. In this micro-canonical simulation 
the temperature and also the specific heat are determined as averages of expressions easy to 
implement. The XK-chain is studied for a test. The second order transition on a cubic lattice 
and the first order transition on fee lattice are analyzed in greater detail to have a more severe 
test about the feasibility of this micro-canonical method. 

PACS numbers: 05.50.+q, 75.10.Hk, 75.40.Mg 

I. INTRODUCTION 

Dynamical simulations of large systems of spins or particles should in principle lead to the same information 
as statistical simulations generating a canonical ensemble. Although the energy is conserved one can determine 
the temperature for a system consisting of kinetic and potential energy the temperature by the average kinetic 
energy. With the recent advance in the understanding of the micro-canonical temperature initiated by Rugh , 
one can "measure" the temperature and also the specific heat along the trajectory of spin systems where such a 
decomposition of the energy does not exist. In the case of XF-spins one could also add the sum over the angular 
velocities but a dynamical check of statistical mechanics by a common molecular dynamics simulation is not an 
easy numerical task . We think by using spin systems alone one has more efficient means to pursue this program. 

Several attempts have been made to establish a strictly micro-canonical approach to study phase transitions |^J^ . 
Here we want to present a study of the planar or Xy-spin systems, making use of a micro-canonical temperature 
definition discussed recently and the "over-relaxation algorithm" ||^,^,|ll|. This algorithm has been used to 
accelerate Monte Carlo simulations for Heisenberg or XY spins for example. Since it conserves energy it can be 
taken as a Metropolis step, where steps leading to equal or lower energy are always accepted. The Metropolis 
procedure leads to thermalisation, whereas the over-relaxation algorithm generates different configurations and 
thereby improves the convergence. The question we want to address is, can one work alone with the over-relaxation 
algorithm. Since the temperature can be calculated as an average the Metropolis procedure is not needed for the 
temperature determination. However, it is not guaranteed whether in a fixed energy simulation one "visits" enough 
regions of phase space in order to determine the physical observables reliably. More specifically we want to ask: 
is the finite size analysis the same as for canonical simulations? A closely connected question: is the analysis of a 
first order transition really different from the standard one using a histogram |^ analysis? 

For numerical tests we have chosen XK-spins on a chain and on three dimensional simple cubic (sc) and face 
centered cubic (fee) lattices both with nearest neighbour coupling. A second order transition on the sc lattice and 
the first order one on fee will thus be studied. 

II. MICRO-CANONICAL TEMPERATURE FOR CLASSICAL AK-SPIN SYSTEMS 

A. General approach 

Instead of three components for the classical Heisenberg system in our preliminary paper here we want to 
^ . study the simpler A"F-spin system where only two components S = (5^, 5^) have to be taken into account, that 
■ " " ' is the spin vector is confined to a circle with radius S. No dynamical equation can be written down similar to 
dS/dt = S X H ior the three dimensional vector 5 in a magnetic field H. The energy —H ■ S = —HScosf is 

the same for the two directions ±ip the vector S can point. Contrary to the three dimensional case, there is no 
continuous path between these equal energy states at ±ip. 

Still we can calculate the size of the phase space W summing over all contributions with the same energy E. W 
corresponds to the density of states in the classical limit. With the entropy 5 as a function of energy given by 
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5 = 111^" (1) 

(Boltzmann's constant set to unity), the density of states contains all information necessary for a thermodynamic 
analysis. We want to show that the derivative of S with respect to the energy is easily accessible. 
With TV spins of length \Si\ — 1, that is Si = (cosfi, sin ipi) one has 

JsiE-nml.d^^^ I ^ (2) 

■H=E 

where the integration is over N unit circles. However, the integral is restricted to 2N—1 dimensions by Ti — E. 
This restriction taken into account, one arrives at the second form of the integral with the gradient taken with 
respect to all ipi angles . Because the direction of the infinitesimal surface elements dO normal to the surface 
of constant energy coincides with the direction of the gradient V7i, one arrives at a second form of the constant 
energy integral 

H=E 

Applying Gauss' theorem this integral can be written as an integral over all phase space with energy lower than E, 
the restriction given by the Heavyside step function 0. The quantity to be integrated is the divergence of a vector 
field X with the property that the scalar product of S/H and X = XlJij jV^p is unity ||^. The derivative is then 

a constant energy integral like (H) for W . The ratio of these two integrals gives then the inverse temperature 

1 dW IdE 
T = — 

since 1/T — dS/dE with the entropy S given by ([l]). 

In order to obtain the inverse temperature according to and (^) one has to average the following quantity, 
called l/Txy for short 

where we use the notation di — d/dipi. To compare with the Heisenberg case j^] of a three dimensional vector 
confined to a sphere one identifies d/dfi as the third component of the angular momentum operator 

— - S'— - S^— (7) 

used in quantum mechanics, only "z" is missing. Eq.^ can then be viewed as a reduced form of a micro-canonical 
temperature for a Heisenberg spin system with the complete angular momentum operator Ci = Si x Vi 

' ^At-^^ (8) 



instead of only the third component of Ci appearing in eq.(^). 

For large system sizes N ^ I the quantity to be averaged in (^) can be simplified. The two terms generated by 
differentiation are 

where the last one can be neglected since it is only of order The reason is the short range of exchange 

interactions. The Hamiltonian Tl connects a spin Si only to its neighbors so that the differentiation with respect 
to ipi is nonzero for I = i and / close to i. As a result the second term scales as oc since a single sum in 
the numerator is outweighed by a double sum in the denominator. Also for large spin systems neglecting again 
1 /A^-corrections one can take the inverse of (0) to determine the temperature T directly as an average of 



J-xy ~ 

that is T = {Txy 



^ Eidfn ' ^^^^ 



B. Micro-canonical temperature in simulations 



Although there exist no dynamical equations which conserve the energy for XY-spin systems, there is a so- 
called "over-relaxation algorithm" in use to generate spin configurations of equal energy. This algorithm is used in 
conjunction with the standard Metropolis procedure, see for example to diminish the effects of critical slowing 
down The Hamiltonian we want to study is 

T-Cex = -^Jij COs{ipi - ipj) , (11) 
ij 

where the sum is over the exchange constants Jy connecting close spins at sites i and j. In an equivalent notation 
this is Hex — — -^ij ' with Si — (cosipi, simpi). The energy Ei contributed by an individual spin Si is 



E, 



with Hex = 5 J2i ■ S^^ total energy Hex one has to sum over all sites of the lattice and since the exchange 
interaction appears twice the sum has to be divided by two. The "local" energy Ei depends on the "molecular" 
field Hi generated by the neighbors. It does not change if 

S^ =^ -S^ + 2H^{S^■H,)/Hf . (13) 

This is the basic step of the "over-relaxation" algorithm which in geometrical terms means a change of Si to its 
mirror position with the mirror line given by Hi as in Fig. 1. In the simulation the spins to be updated are mostly 
chosen randomly ||lj] . 

Using the special form for the Hamiltonian (^) the formula for the temperature (^) simplifies to 



^2 



Txy = - ^'l:^ , (14) 

^ liex 

since summing the second derivative of the Hamiltonian Hex with respect to all angles ipi one recovers Hex merely 
multiplied by —2. The easiest way to test the validity of the temperature determination by averaging (^^ is to 
apply the over-relaxation algorithm to a spin chain. However, this algorithm for a linear chain with 

Hehain ^ ^ J J2iLi cos{(pi - fi+i ) and ipN+i = fi (15) 

would be non ergodic in the case that all exchange constants are equal In alternating the exchange by ±10% 
along the chain one has a remedy for this defect easy to implement. The temperature according to (llj) for a 
homogeneous chain is 

Tchain = -^Z^ill [jJ2±^^'^iVi - V^i±l)]'^ /Hchatn ■ (16) 

The task is then to check whether for a large chain the canonical relation |ic| ] between energy E and temperature 

e^E/N ^~Jh{J/T)/lo{J/T) (17) 

is reproduced, /q and Ji are modified Bessel function. The specific heat as the derivative C — dE/dT then is 

C/N = {J^ +Te-e^)/T^ . (18) 

We checked that the slight change of the exchange constant J = 1.0 ± 0.1 had no visible effect in Fig. 2. For a 
comparison using the average value J = 1 is sufficient. The temperature and the specific heat for a XK-chain 
determined by the over-relaxation method by varying the energy agree well, as can be seen in Fig. 2, where energy 
and specific heat are plotted in the conventional way as functions of the temperature. The question how to calculate 
the specific micro-canonically for XY-spin system will be addressed in the next section. 

An initial state with a prescribed energy to start the simulation is found in the following way. One solves 
numerically the differential equations 

^,=±d,Hex/T. (19) 

A change of the angles (pi of the spins involves of a change of the total energy, that is Hex — ^{^/t) T^ex)^ ■ 

Depending on the size and sign of the relaxation time r in principle starting from randomly chosen orientations 
LPi of spins any prescribed energy could be reached. Very low energies states, however, can only be found by a 
succession of energy conserving over-relaxation steps and steps lowering the energy by small amounts using (19[). 



III. SPECIFIC HEAT OF A XY SPIN SYSTEM 



A. Specific heat for the XF— chain 



In taking the derivative of the entropy S with respect to the energy E one obtains the inverse temperature. The 
second derivative with the respect to the energy will then determine the inverse of the the specific heat. Formally 
the first derivative defines the inverse temperature by the divergence of a vector X 



N 

X = (diH, d2H, ■ ■ ■ BnU) I {d^nf (20) 



according to (||) written as 1/T = diXi. This standard form of X can be changed. One can modify X as long 
as the scalar product between X and V7Y does not change, that is 

N 

Y,X,d,n = 1 . (21) 

i=l 

To get the volume of phase space with (||) correctly, (^ij) is the only condition 1^ . In a "hydrodynamic" picture 
this means that the flow in phase space X across Ti. — E should be the same as the "natural" flow given by ( pO| ) 
which is perpendicular to the TL — E surface. The constraint (|l|) leaves enough room for choices better fitted to 
specific heat calculations as will be shown. 
For a chain a variant to X is 

N/2 2 

X' = (o, d^n, 0, • ■ • 0, d^n) / ^ (d^^n) , (22) 

i=l 

and one easily checks that the normalisation (pi|) is still valid. In the following it is assumed that Ti is of the 
exchange type restricted to nearest neighbors like (|ll|) or specifically (|l^) for the one dimensional case. The inverse 
temperature is then given by an expression similar to 

1 ^ d..n ^ T^&i{H 2^3^11!^ (23) 

where in the summations only the even sites along the chain appear. The second term, although oc l/7V, contributes 
to the specific heat as will be seen. Its simpler form compared to the one derivable from the original expression (^) 
for 1/T is helpful for an easier numerical calculation. 

For the temperature determination the inverse of ( p^ ) without the 1/A^-correction is sufficient, that is we take 

N/2 2 

^ = -E(^2,^) /n (24) 

and determine the temperature by T (T). The energy Ti or E in the denominator is obtained by summing up 
the second derivatives, that is J^iHi '^2i^ = ■ The derivative of T with respect to temperature consists of two 
contributions 

dT _ 1 _ 
dE ~ C ' 

^ /1\ l\^^/'^(F)^.T-l\'^P)2q^\ 

(25) 

The first has its origin in the direct energy dependence of the temperature average of ^ ^ di {X[T) . The second 
term comes from the energy dependence of the normalization necessary for the temperature average. The definition 
( p4[ ) leads to the first two terms in the expanded form of the next line whereas the two last terms are generated by 
using (H) for l/T'. With 

1 ^ 1 ^ 1 _ ^ , ^ .26^ 

T T + AT T ^ ' 




A 



and AT = T—T the '1' in ( |25| ) disappears and only terms cx remain. This is to be expected since the reciprocal 
of the specific heat behaves as dT/dE cx One notices that the inverse of the specific heat depends on the 

temperature fluctuation. As a further simplification the averages of the last term of ( p6| ) can be taken separately 
for the denominator and the numerator by neglecting 1 /N"^ corrections. Further the average of the denominator is 

simply {ETY since (J^ji^^jTi-)^) = T E according to (p^. The difference to the true average would result again 
in a correction and can therefore safely be neglected. 

Taking the simplifications into account, the specific heat is then 

1 T ((AT)^) jE^J'{d2.nYdliH) 

C - E ^ TE^ ^ ^ 

No further simplifications are possible. In the last term appears the local energy — besides the square of the 
"torque" {d2iH)^ used for the temperature determination. To get an understanding of the average of these two 
quantities let us define a local energy deviation /S.Tt2i = —d2{H — 2^TL/N. The 7 value is determined by splitting 
the average of into two parts 

Y.';il{d2,nfdln) = -{Y.^il{d,,Hfi^^^^^ = ^-^te-. (28) 

such that the first term depending on the local energy fiuctuations should vanish. The average of the second term 
is ozT E^ according to (^). The inverse of the specific heat (j2^) is then 

C ^ E^ N T2 
which at low temperatures with T/{~E) <i; 1 should reduce to 

^ ^ £ _ {{ATf) 

c ^ N r2 

a form derivable (see appendix) for a chain of harmonic oscillators with the temperature fluctuation 2/N so that 
Cose — N/2. At low temperatures one expects 7=1, since the XF-spin chain should have the same specific heat 
as a harmonic oscillator chain. 

To get the specific heat of the ATF-chain shown in Fig. 2 relation (^^ has been used. At low temperatures the 
specific heat for one spin is 1/2 as expected. 



(30) 



B. Specific heat for tiie XY— lattice systems 

The strategy to calculate the specific heat for a chain can also be used for two or three dimensional lattices. For 
square and simple cubic lattices a decomposition in two sublattices analogously to the even odd decomposition for 
the chain can be made. A spin is even or odd if the sum of its coordinates (i, j, k) is even or odd. 

The point is that the no nearest neighbor "bond" Si-Sj appears twice in different components which is the case 
for the vector X' defined by (^2|). As a consequence only local energy and temperature terms have to be averaged. 
For a chain and a simple cubic lattice starting from all "even" or all "odd" sites all the interactions Si ■ Sj are 

taken into account and 'Yl!i=i = ^r I = even, odd. So after renumbering the spins ( p7|) can be used 

without further changes. As an example we show in Fig. 3 the specific heat of the ATy-model on a simple cubic 
lattice obtained this way. Only close to the critical region differences between the canonical simulation using the 
Metropolis method and the micro-canonical simulation are visible. 

In the general case, for a face centered cubic lattice as an example, slight changes have to be made. The fee 
lattices can be thought to consist of four sublattices, one simple cubic lattice and three additional ones shifted in 
the three direction (1, 1, 0), (1, 0, 1) and (0, 1, 1) to the centers of the faces of a cube. If one picks out one of 
these simple cubic lattices one can take over all the equations from the last section for the chain, only the energy 
connected to that sublattice is no longer a constant since is not the total energy. With of lattice points of the 
simple cubic lattice, that is TV = 4 the total number the spins, the temperature analogous to ( |2^ ) is 

% = -Y,[d^,in) /Ui with T={Ti), (31) 

i=l 

where the summation takes into account all angles ipj^i of one of the cubic sublattices with number ? = 1, 2, 3 or 4. 
The partial energy is Hi = — X^i^i (^u'^ ■ ^ sense one has four thermometers which should of course give the 
same temperature. Similarly one can calculate the specific heat starting from one sublattice. We studied the 
anti- ferromagnetic transition on a fee lattice (see section V) , but the specific heat calculation with formulas similar 
to ( p7|) did not give stable results in the neighborhood of the first order transition because of the singular nature 
of the specific heat. 



IV. SCALING FOR THE MICRO-CANONICAL TEMPERATURE 



In principle the scaling behavior close to a second order phase transition should be the same for a canonical 
ensemble and for a micro-canonical ensemble. So one expects for a micro-canonical ensemble close to the critical 
energy that the "measured" temperature should dependent on the size of the system. This corresponds to the size 
dependence of the energy for a canonical ensemble at the critical temperature |^ . 

The correlation length ^ grows like ^ cx £,o/\T — T^]'^ for a temperature close to the critical temperature T^. This 
length ^ cannot become larger than the size L of a finite system. In a micro-canonical simulation at the critical 
energy Ec this finite size effect should be directly visible. Since the system size takes the role of the correlation 
length, the recorded temperature should deviate from its critical value Tc according to 

Tl =r, + (32) 

that is by rewriting \Tl ~ Tc\ oc L^^/'^ as an equation. In Fig. 4 the size dependence of the critical temperature is 
plotted and one sees that the temperature defined as an average of ^ reproduces the expected effect. 

We took as the critical temperature E^ = 0.989 found by Schultka et al. In the temperature scaling analysis 
we should recover the temperature these authors started with that is the value Tc — 2.2017 determined by Janke 
p6| . This is actually the case, since we find — 2.202. 

For a canonical system where the temperature is fixed by its critical value Tc the critical energy Ec has finite 
size corrections. Loosely speaking one has \E — Ec\ oc |T — Td^^" in the neighborhood of the critical point. The 
linear dependence between energy changes E — Ec and temperature changes T — Tc is modified to a power law 
dependence by the specific heat exponent a. With this modification the finite size correction of the critical energy 
follows in close analogy to ( p2| ) 

El^Ec + EiL-^^-"'^/-' (33) 



from \El — Ec\ oc The relation (33) has been used to determine the critical energy For an XY- 

system practically the same exponent for the micro-canonical scaling according to (|3^) or the canonical one using 
(33) can be taken, since a is very small Actually in ( |3^ ) a correction term E2 L^^f^^")/'^ has been added to 
(B3). In the analysis of ( ^2| ) the corresponding correction is also needed, as can be seen in Fig. 4. 

We test the over- relaxation algorithm by studying the dynamical scaling behavior To use this procedure 
one should start with a configuration as far as possible from equilibrium with an energy fixed to its critical value. 
We take a random configuration for the spins as a starting configuration. It turns out to have an energy close to 
zero. By moving the spins step by step in integrating the differential equation (^9|) the energy is increased until the 
critical value is reached. Then with the of over-relaxation dynamics the energy stays always at the critical value, 
but the correlation length should start growing like ^ oc t-^^^ as a function of the number of update steps t. Finally 
after many steps the equilibrium situation is reached, that is for ^ « i. Therefore the temperature as a function 



of a moderate number of updates t should follow a power law relation like ( |32D 

Tt^TL + Tit-^/^'-""' . (34) 

using \Tt-TL\ oc 

The magnetization \M\ exhibits more clearly the power law dependence, since no constant term must be sub- 
tracted. It should depend on the update steps t like 

(|M|) oc tT/^^^"^) . (35) 

This corresponds to (M^) oc U^^" or (|M|) oc L^^"^^, valid for finite size scaling with the exponent 7 for the 
susceptibility. The rule of transcription is simply that one should replace L by t^/'' . Since rj is small we take 1 for 
the ratio 7/(2 v) ^ I - rijl. 

We used two slightly different ways for the implementation of the over-relaxation algorithm. The first one 
consists of choosing randomly the spins to be updated. The second method uses a sequential update. The lattice 
is partitioned into two sublattices corresponding to the up and down sublattices of an anti-ferromagnet. Since the 
interaction is only between nearest neighbors the numbering of the spin used for the sequential update does not 
matter. First the the spin vectors sitting on the up-spin sites are updated and then the others. Technically one 
decides whether an index is even or odd before one starts with the over-relaxation. A spin is an even or odd one 
if the sum of its coordinates (i, j, fc) is even or odd. 

Taking the results of the simulation depicted in Fig. 5 one notices that with Zgcq — 1/0.58 « 7/4 the sequential 
update is the more efficient algorithm. However, to be more on the safe side concerning ergodicity one might prefer 
the random update with z^an — 1/0 A2 « 7/3. Since always half the total number of spins are "turned around" the 
greater efficiency of the sequential update might be connected to the nonlocal nature of the update reminiscent of 
the very efficient cluster algorithms. 



V. FIRST-ORDER PHASE TRANSITIONS 



The micro-canonical technique can be shown to be most suitable in simulating first order transitions. A XY-spin 
model having such a transition is an anti-ferromagnet on the fcc-lattice. For simplicity we take only the interaction 
between nearest neighbors. This model has been analysed by Diep and Kawamura using standard Monte Carlo 
techniques. There the first order nature of the transition for XY and Heisenberg spins was shown. Especially for 
an Ising model the frustrated anti- ferromagnetic order on a fcc-lattice leads to a strong first order transition p[p^ . 
For the vector spins this phenomenon is less pronounced. 

The temperature for a large range of energies has been determined and the region between 1.35 and 1.55 was 
scanned more carefully. The energies are all positive and therefore all temperatures negative since we kept the 
convention of the preceding sections that the ferromagnet has the lowest energy and the exchange constant J = — 1. 
In the inset of Fig. 6 the S~shape of the energy versus temperature dependence clearly indicates the first order 
nature of the transition. 

We want to analyse in greater detail the data for the largest lattice with 4 x 24^ spins. The idea is to use 
Maxwell's equal area construction to determine the transition temperature. A direct approach is not advisable 
since the statistical errors for the temperature changes are quite large. A better way is to determine first the entropy 
With the energy density e = E/N, the entropy density s = S/N, the reciprocal temperature (3 = 1/T and 
(3 = ds/de one obtains the entropy by integrating 



The logarithm of the density of states W as the entropy is changed by subtracting E/T, that is P = • e~^^^ . 
This probability P depicted in Fig. 7 has the double peak form of the histograms one obtains by counting the 
frequency an energy E occurs in a Monte Carlo simulation. The peaks must have equal height at the transition 
temperature corresponding to Maxwell's equal area rule. 

The temperature we find by this procedure is Tc = 0.7397. However, in looking at the insert of Fig. a considerable 
one can observe a consiberable size dependence. The energy difference between the two maxima Ae = 0.13 gives 
an estimate for the latent heat. One can also read off from Fig. 7 the difference in the logarithm of the probabilities 
\nPmax — InPmiTi ^ N * 0.0005, that is for InP at the peaks and for the minimum between the peaks. Since 
the cluster has = 4 * 16"^ spins the ratio Pmin/Pmax ~ or 10"'' is already too small for a standard Monte 
Carlo simulation. In upgrading the weights by using the "multi-canonical" method this low probability problem 
is circumvented. With the micro-canonical procedure no modification seems to be necessary here. However, first 
order transitions are beset of hysteresis phenomena, they are bound to occur also in micro-canonical simulations. 



We confirmed that the temperature calculation as averages given by (|l0|) and give reliable results. The severe 
tests were (a) the size dependence of the temperature for a second order transitions and (b) the non monotonic 
dependence of the temperature on the energy in the neighbour-hood of a first order transition. Further we wanted to 
demonstrate that efficient micro-canonical simulations for vector spins are feasible making use of the over relaxation 
algorithm. However, since the histogram technique is missing, the critical energy cannot be determined precisely 
by Binder's cumulant technique p^ . For first order transitions the method appears suitable. 

The inverse of the specific heat as the derivative of the temperature with respect to energy is linked to averages 
of fluctuations. They are far more complicated than the energy fluctuation one has to calculate for the specific heat 
in the canonical ensemble. Nevertheless we could show that the specific heat can be calculated micro-canonically. 
That the temperature definition is not unique is the essential point of Rugh's micro-canonical temperature definition 
P]. This freedom of choice was essential for finding manageable expressions. 



We thank W. Janke for pointing to the close similarity between a micro canonical and canonical analysis of first 
order transitions. We are thankful D. Loison for showing us the use of the over-relaxation algorithm. 
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VII. APPENDIX 



A. Micro canonical Specific heat for a harmonic oscillator chain 



The concept used for the micro canonical temperature and the specific heat determination is best be illustrated 
for a set of harmonic oscillators. We take a one dimensional example with the potential energy 

N 

y = i ^(x; - Xi+i)^ and xn+i=xi (37) 

i 

for a chain of coupled oscillators. Following the procedure for the spin chain of section III a vector X is defined by 

N/2 

X = - (o, F2, 0, F^, ■■■ 0, Fn) I J2 ^2. : (38) 

where the expressions F„ — ~dV/dxn mean the forces acting on the n-th particle. Obviously Rugh's condition for 
the "flow" X, namely — FiXi = 1 holds, since only even components appear in the scalar product. The minus 
sign comes from the definition of the forces Fi as negative derivatives. 
The inverse temperature defined as the divergence of X is then 

with d2i — djdxii. Since the second derivative of the chain potential is &^ — 2 one gets finally 

' (40) 



T' vA'/2p2 



2j 

that is the expression one has to average to obtain the inverse temperature. 

Neglecting the finite size corrections we determine the temperature by averaging the quantity 

NI2 

^ = 4e^2v (41) 



The specific heat C is then given by 



Ar/2 

i=l Z^J = 1 ^2j 

-l + (r)(i,) (42) 

Using the expansion ( p^ ) for the reciprocal temperature the last equation for the inverse of the specific heat reduces 
to 

1 = 1- (43) 

in the limit for 3> 1. This is eq.(|30|) used for a comparison with the specific heat formula ( |2^ ) of the XY-apin 
chain. The last term of ( ^ ) is the temperature fluctuation. Its size must be 2/N since the specific heat of N 
oscillators is C = N/2. One could determine this temperature fluctuation for an oscillator chain analytically, since 
V and T are quadratic forms, but this technical detail would be of little interest. 

The last formula has the same structure as the speciflc heat formula of Lebowitz, Percus and Verlet for a 
system of interacting particles. There the total conserved energy is the sum of the kinetic and the potential energy. 
From the fluctuations of the kinetic energy K (proportional to the fluctuations of the temperature) the specific 
heat C can be determined in the micro canonical simulation with the help of 

i(K-(K)r) ^ A _ 1 (44) 

(i^)2 3iV C ' ^ ' 

where the average of the kinetic energy for N particles "measures" the temperature by {K) — 3NT/2. For the 
oscillator chain discussed here only the potential energy was considered (see (A.l)) and the partition of the energy 
is between even and odd sites. 
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FIG. 2. Specific heat and temperature per spin as a function of temperature for XK-chain of 16 spins. The lines represent 
the canonical relations (14) and (15). 
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FIG. 3. Specific heat of a 24^ cubic lattice, circles canonical and triangles micro-canonical simulations. Near Tc ~ 2.20 
canonical (2- 10"" steps) and micro-canonical results (2-10^ steps) differ. 
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FIG. 4. The temperatures (+) for the critical energy Ec — 0.989 as a function of the lattice size L. One finds from the fit 
(upper curve) = 2.2013 and u — 0.6709 using eq.(|32|) and adding T2 L^^''" necessary to correct the deviation from scaling 
for small L. 
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FIG. 5. \M\ for clusters of size 24^, 32'' and 48^ as a function of sequential updates starting from disorder at the critical 
energy. Also shown is the linear size dependence of |M| used for the scaling plot in the insert. Not shown is the update with 
randomly chosen spins for 32^ with an almost linear relation between \M\/L and (steps)° *'^/L. \M\ has been averaged a 
few thousand times. 
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FIG. 6. Results of microcanonical simulations for fee lattices with linear dimension 8, 12, 16 and 20. The energy is plotted 
as a function of the temperature. The inset shows the S-form of the energy in the vicinity of Tc « 0.75. 
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FIG. 7. Relation between energy density e and probability P{e) = e"^''' for the antiferomagnet on a fee lattice of 

size 16^ aeeording to (Bq). Equal height of the two peaks selects the best critical temperature. 



1 fi 



